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Two centuries of research on phase transitions have repeatedly highlighted the importance of critical fluctua- 
tions that abound in the vicinity of a critical point. They are at the origin of scaling laws obeyed by thermody- 
namic observables close to second-order phase transitions resulting in the concept of universality classes, that 
is of paramount importance for the study of organizational principles of matter. Strikingly, in case such soft 
fluctuations are too abundant they may alter the nature of the phase transition profoundly; the system might 
evade the critical state altogether by undergoing a discontinuous first-order transition into the ordered phase. 
Fluctuation-induced first-order transitions have been discussed broadly and are germane for superconductors, 
liquid crystals, or phase transitions in the early universe, but clear experimental confirmations remain scarce. 
Our results from neutron scattering and thermodynamics on the model Dzyaloshinskii-Moriya (DM) helimagnet 
(HM) MnSi show that such a fluctuation-induced first-order transition is realized between its paramagnetic and 
HM state with remarkable agreement between experiment and a theory put forward by Brazovskii. While our 
study clarifies the nature of the HM phase transition in MnSi that has puzzled scientists for several decades, 
more importantly, our conclusions entirely based on symmetry arguments are also relevant for other DM-HMs 
with only weak cubic magnetic anisotropics. This is in particular noteworthy in light of a wide range of recent 
discoveries that show that DM helimagnetism is at the heart of problems such as topological magnetic order, 
multiferroics, and spintronics. 



I. INTRODUCTION 

Critical phenomena were observed for the first time in 1822 
in the form of the critical opalescence of water vapor when 
Cagniard de la Tour discovered the critical point of the gas to 
liquid phase transition [1 1. The associated change of the phys- 
ical properties is also referred to as continuous (or second- 
order) transition, implying that the order parameter character- 
izing the ordered state, emerges smoothly. Hence, as recog- 
nized by Cagniard de la Tour a critical point represents a "spe- 
cial state" (etat particulier) [2| because the disordered and the 
ordered phase are indistinguishable. 

The unusual physical properties are thereby generally re- 
ferred to as critical phenomena. At the heart of the critical 
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phenomena is an abundance of low-energy fluctuations of the 
order parameter that extend over increasing length scales as 
the critical point is approached. The divergence of their so- 
called correlation length £ results in universal scaling laws for 
observables that only depend on the symmetries of the criti- 
cal system while being independent of its specific microscopic 
details. This led to the notion of universality classes — a cor- 
nerstone of modern physics — providing a common frame- 
work for a wide range of systems with the same critical be- 
havior despite entirely different microscopic character O . 

As one of the most remarkable aspects, it has long been no- 
ticed theoretically that an excess of critical fluctuations may 
change the nature of the phase transition entirely. If the phase 
space available for the critical degrees of freedom is suffi- 
ciently large, the system may evade the critical point to avoid 
the large entropy associated with the fluctuations by realizing 
a discontinuous first-order transition into the ordered phase. 
As a result, the correlation length does not diverge and the 
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order parameter varies also discontinuously at the transition 
which is then accompanied by the release of latent heat. 

Theoretically, fluctuation-induced first-order transitions 
occur as a consequence of non-analytic terms in the 
Ginzburg-Landau free energy functional that are generated by 
fluctuation-induced corrections. Perhaps the best-known ex- 
ample concerns thereby the Coleman-Weinberg effective po- 
tential for an order parameter coupled to a fluctuating gauge 
field 01 . This found important applications in the context of 
phase transitions in the early universe addressing cosmic in- 
flation [5] and the problem of baryogenesis (6). However, 
it is also relevant for the description of superconductors and 
smectic-A liquid crystals (7J. In contrast, for an order pa- 
rameter with a large number of components N > 4, its self- 
interaction may already suffice to drive the transition first- 
order as pointed out by Bak and coworkers |8|. Finally, Bra- 
zovskii [9 1 considered theoretically critical fluctuations that 
become soft not only at a single point in momentum space (as 
e.g. for a ferromagnet), but rather on a finite manifold, notably 
a sphere. In such a case, the density of states for critical fluc- 
tuations exhibits a one-dimensional singularity so that inter- 
action corrections are expected to drive a strong suppression 
of the correlation length and, eventually, a fluctuation-induced 
first-order transition. 

Various settings have been proposed in which Brazovskii- 
type transitions may be expected such as weakly anisotropic 
antiferromagnets [9|, liquid crystals ifTOlfTTI . diblock copoly- 
mers 021, the Rayleigh-Benard convective instability |[T3l . 
pion condensation in nuclear matter lfl4l . and Bose-Einstein 
condensates in multimode cavities ifTSl or with spin-orbit cou- 
pling [16]. However, only few experimental realizations of 
Brazovskii systems such as diblock copolymers ifTTI were re- 
ported so far. The outstanding experimental constraint thereby 
are precision measurements of the relevant fluctuations, which 
are difficult to resolve llT8l . Prior to the present study direct 
experimental evidence for a Brazovskii-transition hence rep- 
resented a major challenge of interest to a wide range of top- 
ics. 

In this paper we identify the onset of helimagnetic (HM) or- 
der in a class of cubic materials as a prime example for a Bra- 
zovskii transition using a quantitative comparison of experi- 
ment and theory. Specifically we consider systems, in which 
the HM order arises from (chiral) Dzyaloshinsky-Moriya 
(DM) spin-orbit interactions in the presence of weak cubic 
magnetic anisotropics. Prominent members of this class of 
DM-HMs crystallize in the non-centrosymmetric P2i3 space 
group, encompassing the metallic B20 compounds MnSi and 
FeGe, the B20 semiconductor Fei-^Co^Si, and the multi- 
ferroic insulator Cu 2 OSe03. These materials have recently 
generated tremendous scientific interest, when a new form of 
magnetic order - a skyrmion lattice - was discovered IF2TI4241 . 
As its novel aspect each skyrmion is characterized by a non- 
zero topological winding number. Moreover, in the skyrmion 
phase of MnSi spin transfer torques have been observed at 
record low current densities j = 10 6 A m~ 2 making these 
materials interesting for spintronic applications ll25l l26l . In 
the insulator Cu 2 OSe03, magnetic skyrmions possess elec- 
tric polarization and the concomitant magnetoelectric cou- 



pling promises interesting multiferroic behavior Il27ll28l . 

In fact, it was argued in Ref. [21] that the skyrmion lattice 
phase only exists due to strong renormalizations attributed to 
thermal fluctuations. This motivates a more detailed study of 
fluctuation effects on the phase diagram in general. In the 
present work, we investigated the influence of fluctuations on 
the HM transition of MnSi at zero field by means of carefully 
designed experiments using a combination of small angle neu- 
tron scattering (SANS), and measurements of the magnetic 
susceptibility and specific heat. Notably, we have mapped out 
the fluctuations in three dimensions demonstrating that they 
indeed emerge on a Brazovskii-sphere. Within Brazovskii 
theory the temperature dependence of the correlation length 
hereby naturally explains all features observed in the thermo- 
dynamic measurements. As our results are entirely based on 
symmetry considerations they are expected to be of general 
importance for DM-HMs. 

Our paper is organized as follows; first, we review the state- 
of-the-art on the HM phase transition in MnSi and the experi- 
mental methods used for this work, followed by a description 
of our experimental and theoretical results, before we con- 
clude with a discussion. 



II. THE HELIMAGNETIC PHASE TRANSITION IN MNSI 

In the absence of a magnetic field MnSi exhibits a HM 
ground state below T c w 29 K ll29l [3D characterized by a 
hierarchy of energy scales due to the weak spin-orbit cou- 
pling Aso- First, the strong ferromagnetic exchange interac- 
tion J aligns the spins on short length scales. Second, the 
lack of inversion symmetry in the P2\i space group allows 
for a weak DM interaction, D ~ 0(Xso)> mat leads to a chi- 
ral twist of the magnetization and stabilizes HM order with a 
pitch of 180 A deep in the ordered phase. Third, weak cubic 
anisotropies that arise in higher order in Aso finally lock the 
magnetic helix in a cubic (111) direction |29l . Importantly, 
this hierarchy of energy/length scales is also reflected in the 
nature of critical fluctuations as T c is approached. For T 3> T c 
the correlation length £ is short and the fluctuations have an 
essential ferromagnetic character. However, as the tempera- 
ture is lowered and £ reaches the order of the DM length scale 
£dm ~ J/D, the fluctuations start to accumulate uniformly 
on a sphere in momentum space of radius Q = D/J and, as 
a consequence, the magnetic correlations develop an oscillat- 
ing character. Finally, as the correlation length increases even 
further, £ > £ cu b, cubic anisotropies favor the fluctuations to 
carry momentum in the crystallographic (111) directions. 

While the HM transition is expected to be second-order on 
a mean-field level, interactions between the HM fluctuations 
were theoretically predicted to give rise to important correc- 
tions driving the transition first-order. The precise mecha- 
nism, however, depends crucially on the strength of the in- 
teraction that generates an additional scale, i.e., the Ginzburg 
length £ci P2l . see Fig.[TJa). According to the Ginzburg cri- 
terion, interactions can be treated perturbatively in the case of 
short correlations, £ <C £gi> while the fluctuations are strongly 
interacting if £ > The limit of very weak interactions, 
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FIG. 1 : Helimagnetic (HM) phase transition of MnSi at T c = 29K investigated by small angle neutron scattering (SANS), (a) There are two 
crossovers for non-interacting fluctuations in a Dzyaloshinskii-Moriya (DM) HM as T c is approached and the correlation length £ increases. 
Interactions result in an additional Ginzburg length £oi an d depending on its size different scenarios are realized (see text), (b) The SANS 
setup used to study the temperature-dependence of £, where the neutron beam scatters off the magnetic helix. Turning the sample around the 
rotation axis parallel to the [112] zone axis (red arrow) the sample can be rocked through the Bragg condition for the magnetic helix. With 
this sample orientation the magnetic Bragg condition can be fulfilled for the two magnetic propagation vectors Qi = [111] and Q2 = [111] 
[white arrows in panel (c)] corresponding to two of the four possible helical Q-domains (see text), (c)-(j) Magnetic intensity distribution of 
MnSi close to T c . Below T c [panels (c)-(e)] discrete magnetic Bragg spots corresponding to the helical order are visible, whereas above T c 
[panels (g)-(j)] the magnetic intensity spreads out over a sphere in reciprocal space. The background scattering was determined well above 
the magnetic phase transition and was subtracted from all data sets. The black broken lines in (c) indicate weak reflections due to multiple 
scattering effects along the [001] and [110] directions. Note, that the intensity in panels (c)-(f) is plotted on a logarithmic scale for a better 
comparison between the intensities on the helical satellites and the sphere. 



£ci 3> Ccub 3> £dm> was considered by Bak and Jensen 
||29l who argued that in the regime £ 3> £ cu b an effective 
field-theoretical description emerges in terms of an order pa- 
rameter with N = 8 components corresponding to amplitude 
and phase of helices with momenta along the four equivalent 
(111) directions. The residual interaction between order pa- 
rameter fluctuations might then drive the transition first-order. 
On the other hand, if the cubic anisotropies are associated 
with the smallest energy scale, i.e., the largest length scale, 



£cub 3> £ci 3> (dMi HM fluctuations are already strongly in- 
teracting while they are still uniformly distributed on a sphere 
in momentum space. This scenario coincides with the situa- 
tion considered by Brazovskii |9|, and an interaction-induced 
first-order transition preempts the cubic crossover in this case. 
For the zero-temperature HM transition that is observed in 
MnSi as function of pressure, the corresponding Brazovskii 
scenario was studied theoretically by Schmalian and Turlakov 
ll30ll . Finally, a third scenario arises for very strong interac- 
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tions, £ cu b 3> £dm ^ £cii where the fluctuations are already 
strongly interacting before they develop a preferred chirality. 
Here, the physics on length scales > £ > £ci is gov- 
erned by the Wilson-Fisher renormalization group fixed-point 
resulting in a universal DM crossover at £ ~ £dm- We demon- 
strate below that for MnSi the Brazovskii scenario is realized. 



Although MnSi has been the topic of intense scientific study 
for decades, the nature of the transition from the paramag- 
netic (PM) to the HM phase at zero field is still debated con- 
troversially. Initially, the transition was interpreted to be of 
second-order based on the specific heat in polycrystals 11331 . 
an abundance of paramagnon fluctuations also observed in 
neutron scattering [34-36], as well as a pronounced Curie- 
Weiss behavior in the magnetic susceptibility up to T c (37). 
However, a careful analysis of the temperature dependence of 
the HM Bragg peaks suggested that the transition has a weak 
first-order character [38]. The discovery of partial magnetic 
order in MnSi at high pressures [39 1 inspired a search for spin 
textures with non-trivial topology. In particular, the observa- 
tion of a broad shoulder at T* « T c + 1 K, revealed in more 
detailed measurements of the specific heat ||40| . led to the the- 
oretical proposal of a skyrmion liquid phase as a generic pre- 
cursor phenomenon at essentially each PM to HM transition 
in DM-HMs fiTI . This interpretation was adopted by Pappas 
et al. ll42l l43l who argued that the so-called chiral fraction 
determined via polarized neutron scattering studies provides 
evidence for a skyrmion-liquid phase. Similarly, Hamann et 
al. qualitatively explained their neutron scattering results as a 
complex form of long-range order, a so-called magnetic blue 
phase, motivated by an amorphous lattice of skyrmions ob- 
served theoretically 1441 . However, on closer inspection all 
claims of a skyrmion liquid are highly debatable since they ei- 
ther involve unconventional terms in the Landau theory, do not 
account for the relationship of chiral fraction with skyrmion 
liquid or, last but not least, were obtained with Monte Carlo 
calculations on relatively small system sizes, respectively. 



As an alternative, several studies suggested that the unusual 
specific heat anomaly in MnSi can be explained in the tradi- 
tional framework of HM fluctuations in the presence of cu- 
bic anisotropics. For instance, in 2007 Stishov et al. |45l 
demonstrated by means of the specific heat, thermal expan- 
sion and electrical resistivity that the HM transition is indeed 
first-order with a tiny latent heat. These authors emphasized 
that the broad shoulder in the specific heat, noticed earlier 
|40l|4l]], might be caused by HM fluctuations. Yet, the ori- 
gin of the first-order transition remained unresolved. Further, 
Bauer et al. |06| showed that the specific heat close to T* 
exhibits a so-called Vollhardt invariance [47 1, i.e., a crossing 
point that is invariant under small magnetic fields. Finally, 
Grigoriev et al. [48] interpreted a corresponding signature at 
T* in the magnetic susceptibility and the correlation length 
as a crossover in the character of the fluctuations as the HM 
transition is approached. 
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FIG. 2: (a)-(c) Dependence of the magnetic fluctuations on the rock- 
ing angle ui up to ±35° are exemplaryily shown for T = 29.6 K. 
The sample was rocked around the [112] -axis [see Fig. |TJa)]. (d) 
Integrated intensities corresponding to the red integration regions in 
panel (a) as a function of uj for different temperatures above and 
below T c . The integration regions were selected such that all in- 
cluded points are rotated through the detector plane with approx- 
imately identical velocity. Further, the integrated intensities have 
been corrected for the change of neutron absorption due to the in- 
creasing path length through the sample for increasing angles. The 
red and green shaded regions denote the integrated intensity for the 
critical magnetic fluctuation at 29 and 29.6 K, respectively. 



III. EXPERIMENTAL METHODS 

The small angle neutron scattering (SANS) experiment was 
performed using neutrons with wavelength A = 9.7 A on 
the beamline MIRA [4-9 1 situated at the Forschungsneutronen- 
quelle Heinz Maier-Leibnitz (FRM II). The instrumental res- 
olution was selected by a computer-controlled variable source 
aperture of rectangular cross-section that was installed after 
the monochromator and a cadmium aperture of approximately 
15 mm diameter in front of the sample and is comparable with 
that of other studies ||2T1|431 . The direct beam was masked by 
a Cd mask in front of the detector. For the experiment the sin- 
gle crystal was cooled with a closed cycle cryostat with sam- 
ple tube (CCR) available at FRM II. The single crystal used 
for the SANS measurement is a disc of diameter ps 20 mm 
and 2 mm thickness and was cut from a ingot that had already 
been studied before in several neutron experiments and was 
grown by the Bridgeman technique. l35l l50ll . The residual 
resistivity ratio (RRR) of this sample is m 80. 

The specific heat and ac susceptibility were measured with 
a Quantum Design Physical Property Measurement System at 
temperatures down to 2 K and in magnetic fields up to 9 T. The 
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FIG. 3: Radial q-scans through the magnetic fluctuations above T c = 29 K are shown for all measured temperatures. The inset in each panel 
shows the corresponding direction of the scan in reciprocal space. The q-scans were extracted from the two-dimensional detector images in 
Fig.[T|by performing radial bins with an azimuthal width corresponding to the experimental resolution. The solid lines are fits to Eq. {T} (see 
text). The grey shaded region shown in each panel denotes the experimental resolution. 



specific heat was measured with a standard heat-pulse method, 
where heat pulses were limited to ~ 0.2 % to prevent thermal 
smearing of the signature of the first order transition. For both 
measurements the magnetic field H was applied along a (1 10) 
axis. The single crystal employed for the bulk measurements 
was grown by optical float zoning as detailed in Ref. Il46ll and 
was cut into a cuboid shape (2x2x1.5 mm 3 ) with the two 
large surfaces aligned perpendicular to a (110) axis. The cor- 
responding mass of this sample is 36.75 mg and the RRR is 
80. 



IV. SANS CLOSE TO THE HELIMAGNET TRANSITION 
OF MNSI 

We start with a qualitative discussion of our SANS mea- 
surements before turning to a quantitative analysis. Within 
the ordered phase, the helices carry momentum Q parallel to 
one of the four cubic (111) directions resulting in the forma- 
tion of four domains. In our experiment [Fig.[TJb)] the sample 
was oriented such that the Bragg condition is fulfilled for the 
two pairs of HM Bragg satellites associated with the domains 
Qi || [111] andQ 2 || [111] [white arrows in Fig. [TJc)]. Below 
T c only the discrete Bragg satellites are visible [Fig.[TJc)-(e)], 
but when the temperature approaches T c [Fig. [TJe)] the mag- 
netic intensity drops significantly and the Bragg spots start to 
broaden azimuthally. For T >T C magnetic intensity emerges 
on a ring with a radius that corresponds to the modulus, Q, of 
the helical propagation vector [Fig. |TJf)-(j)] . 

To verify the Brazovskii scenario, we have determined the 
magnetic intensity distribution in three dimensions of momen- 
tum space via rocking scans that were carried out by rotating 
the sample around its vertical axis [Fig.[TJb)]. In Fig.|2|a)-(c), 
we show representative detector images for rocking angles uj 
up to 35° for T =29.6 K. The magnetic intensities for all mea- 
sured rocking angles are similar, showing that the fluctuations 
indeed emerge on a sphere. The magnetic anisotropy on the 



sphere is investigated in Fig.[2|d) which shows integrated in- 
tensities as function of ui for several temperatures below and 
above T c . The intensity increases significantly at large ui as 
soon as the critical magnetic fluctuations arise at T c [Fig.[2jd), 
red triangles]. Here the sharp peak at the center is due to re- 
manent HM order (red triangles, see discussion below), but 
also the intensity of the fluctuating part [Fig.[2|d), red shaded 
region] initially exhibits a shallow maxima at the position of 
Qi (to = 0) due the cubic magnetic anisotropy. However, as 
T is further increased the distribution of magnetic intensity on 
the sphere becomes more and more isotropic [Fig. [2jd), green 
shaded region]. 

In Fig.Qa) we show the integrated intensity of all four do- 
mains (black squares) that measures the order parameter of the 
HM phase as a function of temperature (see SI Text on how the 
intensities were integrated). The magnetic intensity features a 
sharp drop of about two orders of magnitude at T c , showing 
that the phase transition is indeed first-order. For compari- 
son we have also plotted the intensity of the critical fluctua- 
tions just above the transition (blue circles) demonstrating that 
the magnetic intensity of the Bragg peaks at low temperature 
measuring the static HM order parameter is entirely recov- 
ered above T c in the form of critical fluctuations uniformly 
distributed on a sphere in momentum space. 

For a quantitative analysis of the SANS data we have per- 
formed fits of the radial q-scans shown in Fig. [3] Importantly, 
we find that the magnetic intensity at T > T c appearing on the 
sphere in momentum space can be quantitatively explained in 
terms of critical HM fluctuations. The radial q-scans are well 
accounted for [see Fig. |3ja)-(d)] by the modified Lorentzian 
profile as derived by Grigoriev et al. ll36l 
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FIG. 4: Characteristics of the Brazovskii-type fluctuation-induced first-order transition of MnSi determined by SANS. Parameters in panels 
(b)-(d) were obtained from the magnetic intensity by means of Eq. l|T]) (cf. Fig. [3j. (a) The magnetic intensity on the helical Bragg reflections 
(black squares) exhibits a discontinuous jump at T c as expected for a first-order transition. Above T c the magnetic intensity of the helical 
state is recovered entirely in form of critical magnetic fluctuations arising on a sphere in momentum space (blue circles), (b) Magnitude of 
the helical propagation vector Q in the HM phase (black squares) and the radius of the sphere above T c (blue circles) vs. temperature T. (c) 
Cubic anisotropy aj? ub vs. T. (d) The inverse magnetic correlation length k — 2ir/£ vs. T (blue circles). The red markers denote k according 
to previously published results 1361 . The blue solid line is a fit to Eq. {3} that describes the renormalization of k(T) for T — > T c as expected 
for a Brazovskii transition (see text). The black solid line highlights that for T 3> Tmf mean-field behavior is obtained. The black circles are 
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separate regimes above T c ; (i) for £ < £dm [k > Q, see (c)] fluctuations are essentially ferromagnetic as previously shown in Refs. 1341 1351 ; 
(ii) they develop an isotropic chiral character for £ > £dm as shown in the inset of (d); (iii) for £ > ^Gi [k < «Gi] the interaction suppresses 
the transition temperature resulting in a fluctuation-disordered regime just above the fluctuation-induced first-order transition at T c . We find 
for MnSi T c « 29 K and T MF « 30.5 K and thus a suppression of AT = T MF - T c w 1.5 K. 



TABLE I: Summary of the fit parameters obtained by fitting the in- 
verse correlation length k(T) determined via SANS [see Fig. |4jc)] 
and the magnetic susceptibility x(T) [see Fig. [5|a)] by means of 
Eq. (|3j and Eq. (H), respectively. 
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Here q is the reduced wave vector with q = |q| measured 
from the neighboring reciprocal lattice vector, a cu b measures 
the cubic anisotropy, Ub is the Boltzmann constant, and A 



is a proportionality factor that depends e.g. on the magnetic 
form factor of the Mn ions. Eq. ([TJ describes the intensity of 
critical magnetic fluctuations with inverse correlation length 
k = 27r/£ emerging on a sphere with radius Q in reciprocal 
space. However, depending on the magnitude of the magnetic 
anisotropy a cu b the sphere will have shallow maxima along 
the (111) directions due to the cubic invariant (q* + (iy + qi~ 
1/3), where q — q/q. 

We employed simultaneous fits of q-scans along the [111], 
[111], [001] and [110] directions to determine the temperature 
dependence of a cu b- This is in contrast to Ref. [36] where 
a cu b was fixed, resulting in an inverse correlation length k de- 
pending on the direction close to T c in contradiction to Eq. ([TJ) 
[cf. red crosses and triangles in Fig. |4fd)]. Our fits describe 
the intensity observed on the entire sphere (not only along 
the fitted directions) remarkably well for all observed tem- 
peratures as shown in the SI Text. For T = T c + 0.1 K and 
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T c + 0.2 K the fits improve significantly by adding a Gaussian 
profile to Eq. ([T} for the scans parallel to [111] and [111] indi- 
cating the presence of a tiny fraction of HM order above T c . 
This suggests that remanent droplets of the HM phase survive 
above T c as expected for a first-order transition lfl8l . 

The results of the fits are summarized in Figs. |4|b)-(d). Q 
does not show any pronounced anomaly close to T c and keeps 
increasing as function of temperature with about the same rate 
as observed in the HM phase [Fig.Qb)]. In the HM phase Q 
was determined via fits of the helical satellites using a Gaus- 
sian profile (see SI Text). At T c , we find Q = 0.039 A" 1 , 
yielding a pitch length £dm = 2n/Q = 160 A. 

The cubic anisotropy a cu b is shown in Fig. |4|c); whereas 
it is negligible at higher temperatures, it reaches a maximum 
value of a^ uh « 0.023 at T c . This gives a cubic length scale 
£cub = 27r/(a cu bQ) ~ 1040 A. As explained above, the two 
lengths £dm and £ cu b give rise to crossovers in the character 
of the critical fluctuations, see Fig. |TJa). Whether the latter 
cubic crossover develops or not depends on the size of the 
correlation length at the first-order transition. 

The inverse correlation length k is shown in Fig. |4fd) and 
is well in agreement with previously reported results [36]. 
It assumes a finite value at the transition that is given by 
£ c = 2tt/k c s» 1260 A so that it is in fact of similar magni- 
tude as the cubic length, £ cu b- This implies that the first-order 
transition takes place before the cubic anisotropies have fully 
developed clearly disfavoring a Bak-Jensen scenario. 

This is borne out by the value of the Ginzburg length which 
is at the origin of the peculiar temperature dependence of k 
as we demonstrate in the following. The temperature depen- 
dence of k was previously fitted with an ad-hoc two-stage 
power-law dependence l48l . We find, however, that it can be 
naturally explained in terms of a Brazovskii renormalization. 
If the singular fluctuation corrections are taken into account 
self-consistently, one obtains the Brazovskii equation for the 
inverse correlation length [9| (see SI Text for details) 



l MF 



"Gi 



(2) 



where kmf oc T — Tmf measures the distance to the mean- 
field transition temperature Tmf, and KQi can be identified 
with the inverse Ginzburg length, £ci = ^tt/kqi. The explicit 
solution to this cubic equation reads 



k(T) 



/ (t+(1-t3 + v /1-2t3)i/3)2 
KGi V 2V3(i _ r 3 + (Vl - 2r 3 )V3) ' 



where r = (2 1 / 3 /3)k mf /^; = (T - T M f)/T . As shown 
in Fig.|4|d) (blue solid line) the inverse correlation length that 
was experimentally determined from our SANS data in the 
temperature range between T c and 30.5 K is perfectly de- 
scribed by the Brazovskii equation ([3) with fit values listed 
in Table U 

The resulting Ginzburg length amounts to £ci ~ 330 A 
demonstrating that the three length scales obey (dm < £,a < 
( cu b, see Table |lT] and also Fig.|4|e). In turn, the Brazovskii 
scenario does indeed hold. At large temperatures T 3> Tmf, 



i.e., £ <C £ci, the system is disordered already on the mean- 
field level, and Eq. ([3]l recovers the mean-field behavior k « 
%f, as shown by the black solid line in Fig.Qc). However, 
for T ps Tvif, i-e-, £, ~ £gi> the strong one-dimensional sin- 
gularity of the interaction correction prevents the correlation 
length to become infinite and thus impedes the condensation 
of long-range order. Consequently, the transition temperature 
is suppressed by AT = Tmf — T c m 1.5 K, see Table [I] and 
Fig. |4| e ), giving rise to a fluctuation-disordered regime for 
T c < T < Tvif [blue shaded region in Figs. Qd) and (e)]. 
Finally, the fluctuations trigger a first-order transition into the 
HM state at a critical t c < with a numerical value of order 
one, |r c | ~ 0(1), whose description is beyond Eq. |2]i but is 
explained in detail in the SI Text. 

TABLE II: Estimate of various length scales close to the helimag- 
netic transition of MnSi, see text 

chiral DM length £ DM » 160 A 

Ginzburg length £gi ~ 330 A 

cubic anisotropy length £ cu b ~ 1040 A 

correlation length at T c £ c w 1260 A 



V. CRITICAL THERMODYNAMICS OF MNSI 

We now turn to an analysis of the thermodynamic measure- 
ments. The longitudinal magnetic susceptibility for temper- 
atures T > T c and zero magnetic field H = can also be 
expressed in terms of the inverse correlation length 



X\T>T a 



Xo 



1 + k(T) 2 /Q 2 ' 



(4) 



where xo is a constant, and Q is the length of the helical prop- 
agation vector as before. The inverse of the measured mag- 
netic susceptibility for H = is shown in Fig. |5ja) together 
with a fit to the Brazovskii formula of Eq. ([3]), which perfectly 
accounts for the observed T-dependence. The fit parameters 
are given in Table [I] and agree well with the values extracted 
from SANS. To highlight the relationship between k and x, 
we have added the latter in Fig.Qd) (black circles). 

The Brazovskii renormalization of k leads to a striking 
modification of the temperature dependence of the suscepti- 
bility. On the mean-field level, kmf oc T — Tmf, the suscepti- 
bility shows Curie- Weiss behavior and would keep increasing 
as the temperature is lowered [cf. red solid line in Fig. [5ja)] . 
However, the renormalization of n weakens this increase and, 
in particular, induces a turning point 0^x\t=t* = at a 
temperature T* = 30.8(1) K, that is slightly larger than T M f 
(green solid line). To make contact with previous work |48l, 
we also note that d\\ exhibits a local minimum at a temper- 
ature T min = 29.8(2) K with T* > T min . Within Brazovskii 
theory, both these temperature scales are determined by the 
Ginzburg length £q; (or, more precisely, by the ratio (gi/^dm) 
and thus originate from the interaction between chiral fluctu- 
ations. 
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FIG. 5: Signatures of the Brazovskii scenario in MnSi as reflected in the thermodynamic observables. (a) The inverse ac-susceptibility 
(x — x')~ m zero magnetic field vs. temperature T. Here \ is a small contribution to \ that depends on exact position of the sample in 
the measurement setup. The blue solid line is a fit of the magnetic susceptibility to Eq. ptt, T c — 29 K has been determined as the minimum 
in x _1 - The red solid line shows the (field-cooled) magnetic susceptibility in the mean-field approximation. The green solid line displays 
9tX(T) indicating the position of a turning point in x(T) 0tX(X) = 0, horizontal dashed line) at T* ~ 30.8 K. T m i n « 29.8 K denotes 
the minimum in df>x(T). (b) The magnetic susceptibility x~ x' vs - T measured for magnetic fields H — and 0.05 T, respectively. The 
blue solid line denotes the magnetic susceptibility calculated via self-consistent Hartree-Fock-Brazovskii (HFB) theory (see text), (c) The 
specific heat C vs. T for different magnetic fields H, Here T c and T* indicate the position of the sharp first-order like feature and a Vollhardt 
invariance for which 8hC(T) = 0, respectively. We note that the Vollhardt invariance coincides with the turning point in x- The solid lines 
are guides to the eye. The inset shows the magnetic phase diagram of MnSi as determined via x with H applied along [110]. The green squares 
are the positions of the turning point in x(T) that coincide with the crossing point in C in the limit H — > 0. HM, FP, FD, and PM denote the 
helimagnetic, the field-polarized, the fluctuation-disordered paramagnetic and the mean-field disordered paramagnetic regimes, respectively. 
A labels the A-Phase that shows the skyrmion lattice I; 211 . (d) The magnetic contribution to the specific heat C m a g for three magnetic fields 
H vs. T, obtained by subtracting the electronic and lattice parts (see text). The solid lines represent calculations for C mag according to 
self-consistent HFB theory. Thermodynamics does not show clear signatures at the temperature Tmf itself where £ « ^Gi (cf. Fig. [4]), but 
rather features in x an d C at T* w 30.8 K whose value is determined by the ratio £gi/£dm (see text). 



As explained in the SI Text, the self-consistent Hartree- 
Fock-Brazovskii (HFB) approximation goes beyond Eq. <|3j as 
it also describes the fluctuation-induced first-order transition 
and the resulting behavior within the ordered phase. Although 
a self-consistent Hartree-Fock theory is inadequate for the 
description of second-order transitions \5l\, it is controlled 
for Brazovskii systems in the limit of weak interactions, i.e., 
£DM/£ci < 1. However, as we find £oM/£ci ~ 0.5 for 
MnSi, this condition is not well obeyed so that quantitative 
corrections to HFB theory can in principle be sizable. 

In Fig. |5jb) a comparison between experiment and theory 



is shown for the zero-field cooled (ZFC) magnetic suscepti- 
bility. The magnetic susceptibility within the ordered phase, 
T < T c , is anisotropic and thus depends on the domain popu- 
lation. In case of zero-field cooling where all (111) domains 
are equally populated this anisotropy results in the decrease 
of x for T < T c giving rise to the peculiar hooked shape that 
is nicely reproduced by the HFB approximation. The value 
of T c , however, is slightly overestimated by theory and the 
saturation value of the ZFC susceptibility within the ordered 
phase differs by approximately 10%. A small field tends to 
align Q with the field direction resulting in an enhancement 
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of the susceptibility at H = 0.05T in Fig. 5(b) as compared 
to zero field. 

The specific heat C of MnSi as measured for different mag- 
netic fields H is shown in Fig.[5jc). At zero field a sharp first- 
order spike is observed at T c . With increasing H it shifts to 
slightly lower temperatures, while simultaneously its magni- 
tude is suppressed, well in agreement with the known phase 
diagram shown in the inset of Fig. |5Jc). The first-order peak 
is accompanied by a broad shoulder with a slope that de- 
creases with H resulting in a characteristic crossing point at 
T* = 30.8 K. At large fields H > H c2 0.55 T 06|, the 
HM order is suppressed and the first-order transition disap- 
pears while the shoulder develops into a broad feature indicat- 
ing the crossover from a PM into a field-polarized regime. 

Crossing points of the specific heat in general have been 
discussed by Vollhardt fl7l who pointed out that they are 
linked to inflection points of a certain conjugate variable by 
a Maxwell relation. In our case, the crossing point in C is 
related to an inflection point in the magnetization M, 

= d H C\ T , = Td$M\ T , H £° THd%x\ T »- (5) 

In the limit of small fields where M — \H, this is equivalent 
to a turning point of the susceptibility \ which, as discussed 
above, is induced by the Brazovskii renormalization of the 
correlation length. The turning point in \ an d the crossing 
point in C are thus different manifestations of the same phe- 
nomenon. The strong renormalization effects arising from the 
interaction among chiral fluctuations is therefore also at the 
origin of the Vollhardt invariance observed in C. 

Fig. [5jd) shows the magnetic contribution to the specific 
heat C mag that has been obtained by subtracting the electronic 
and lattice contribution using the values established in Ref. 
Il46l . After having fixed all parameters by a fit to the sus- 
ceptibility, the self-consistent HFB approximation predicts the 
behavior for C mag as shown by the solid lines. It nicely ac- 
counts for the shoulder and the Vollhardt invariance close to 
T* and explains the quasi-saturation just above the first-order 
transition even though the parameter ^dm/Cg^ that controls 
this approximation, is not particularly small. The HFB ap- 
proximation becomes less accurate with increasing field and 
at H = 0.2 T ceases to be reliable (see SI text). Within the 
ordered phase, T < T c , the HFB approximation shows a pro- 
nounced additive contribution to the specific heat reminiscent 
of the jump in C on the level of mean-field theory. Inter- 
estingly, this offset is however not observed experimentally 
which might be attributed either to an additional higher order 
A/ 6 -term in the Ginzburg-Landau expansion or to the influ- 
ence of HM Goldstone modes [ 50 1 (see SI text). 



VI. DISCUSSION AND OUTLOOK 

In summary, we showed that the transition from the PM 
to the HM phase in MnSi is driven weakly first-order by 
fluctuations. The first-order nature of the transition is con- 
firmed by our SANS measurement, and the intensity of the 
critical fluctuations for T > T c is found to be practically 



uniformly distributed on a sphere in momentum space with 
only small corrections due to cubic anisotropies. This identi- 
fies the Brazovskii scenario as the driving mechanism for the 
fluctuation-induced first-order transition and rules out the sce- 
nario proposed by Bak and Jensen. This is corroborated by 
the observed suppression of the transition temperature and, 
in particular, by the temperature dependence of the measured 
correlation length that can be quantitatively described by the 
lowest-order self-consistent Brazovskii approximation. We 
discussed the relation between the critical fluctuations and 
thermodynamic quantities, and, especially, explained how the 
Brazovskii renormalization of the correlation length induces 
an inflection point in the magnetic susceptibility and a Voll- 
hardt crossing point in the specific heat. The provided evi- 
dence for these thermodynamic signatures convincingly rules 
out alternative explanations in terms of intermediate skyrmion 
liquid phases as proposed in Refs. HTH441 . 

Our study also offers interesting perspectives for future re- 
search. An obvious question to be addressed is the evolution 
of the Brazovskii fluctuations with increasing pressure. The 
transition temperature of MnSi can be suppressed by applying 
a pressure of p w 15 kbar accompanied by the emergence 
of partial magnetic order and an unusual extended non-Fermi 
liquid phase |39ll52l . The possible relation of these phenom- 
ena to a quantum version of Brazovskii theory 11301 is an in- 
teresting open question. Furthermore, as the magnetic field 
increases the first-order nature of the transition weakens and 
it eventually becomes second order resulting in a tricritical 
point that is to be investigated [cf. inset of Fig.Bfc)]. 

As our arguments are derived from basic symmetry princi- 
ples we expect that cubic DM-HMs with weak anisotropies 
like the B20 compounds generically exhibit a fluctuation- 
induced first-order transition of Brazvoskii type. However, 
one may raise the question whether such transitions also pre- 
vail in DM-HMs with lower symmetry, such as Ba 2 CuGe 2 07 
(tetragonal) |53| - for which critical magnetic fluctuations 
similar to the ones in MnSi have been already observed l54l . 
Ba3NbFe 3 Si 2 Oi4 (trigonal) El, andNdFe 3 (B0 3 ) 4 (trigonal) 
11531 . Notably, DM-HMs with lower symmetry are relevant for 
multiferroics ||57| and magnetic sur- and interfaces ll58H60l 
that are promising systems for memory/storage devices. 

Our work identifies the cubic DM-HMs as a class of sys- 
tems that realize fluctuation-induced first-order transitions of 
the Brazvoskii type, whose unequivocal experimental confir- 
mation has proven elusive in the past. This paves the way for a 
more detailed experimental study of such transitions and their 
properties, e.g., the nucleation of critical HM droplets |6T1 . 
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SUPPORTING INFORMATION 

In the first section of this supporting information details of 
the analysis of the small-angle neutron scattering (SANS) data 
are presented. The remaining sections explain the theoretical 
analysis. In the second section the Ginzburg-Landau theory 
for cubic helimagnets (HM) is introduced and the mean-field 
properties are reviewed. The third section discusses the appli- 
cation of the Hartree-Fock-Brazovskii (HFB) approximation 
to helimagnets, and the fourth section presents estimates of 
the model parameters and energy scales for MnSi. 



I. ANALYSIS OF THE SANS DATA 
A. Integration of magnetic intensities 

We have separated the intensities belonging to the helical 
Bragg satellites and the magnetic fluctuations observed on the 
sphere by carefully choosing appropriate integration regions 
as shown in Fig. [6] Here the segments with the red solid and 
white broken border were used to integrate the intensities for 
the helical satellites and the ring, respectively. 

The segments for the ring have been broken up into six 
pieces in order to remove the intensity of the peaks that de- 
velop due to double scattering in the sample (see black bro- 
ken line in Fig. [6]). The full integrated intensity of the ob- 
served ring was consequently obtained by assuming that the 
intensity is isotropically distributed on the ring and scaling 
the intensity by it/In, where Ir is the combined arc length of 
all six segments. The intensity on the whole sphere can be es- 
timated from the intensity of the observed ring by considering 
the angular resolution of the experiment Soj = 4.5° = 0.0257T 
(FWHM), that was extracted from Gaussian fits to rocking 
scans over the magnetic satellites (cf. Fig.|2]in main text). We 
note that cold triple axis measurements on this sample show 
that the magnetic satellite reflections are narrow (< 0.1°) 0J, 
thus demonstrating that the measured peak width in our ex- 
periment is resolution limited, and the angular resolution in 
a rocking scan is characterized by Soj. The surface segment 
of the sphere covered by the ring at a single position of oj is 

S = 2 Jq sm(8)d6 dip — 4Soj, indicating that the inte- 
grated intensity of the sphere should approximately amount to 

1 Sphere = ^ T R = 40Ir, where I R is the intensity of the 
ring. I sphere is plotted as the blue circles in Fig. |4|a) in the 
main text. 

1 2 

The integrated intensity on the helical satellites I s ' at be- 
longing to domains Qi and Q2 has been calculated by sub- 
tracting the contribution from the ring to the intensities in the 
red solid annular segments (I RS ): I Sat = I RS - IrIrs/k, 
where Irs is the combined arc length of all four red segments. 
The integrated intensity for all four helical domains 7|" t has 
been consequently obtained by multiplying Ig t by a factor 

2 based on the knowledge that all four domains are generally 
equally populated (T). Ig£ t is plotted as the black squares in 
Fig.Qa) in the main text. 



B. Fits of the helimagnetic Bragg peaks below T c 

To obtain the temperature-dependence of the magnitude of 
the helimagnetic propagation vector Q we have performed fits 
of the radial q-scans through the helimagnetic Bragg peaks 
shown in Fig. [7] The peaks were fitted with a Gaussian line- 
shape where the peak position determines the distance to the 
center of the Brillouin zone and therefore the magnitude of Q. 
The temperature dependence of Q is shown in Fig.Qb) in the 
main text. 



C. Calculation of the intensity of magnetic fluctuations 

The fits of the magnetic fluctuations above T c to Eq. ([TJ in 
the main text have been carried out by doing combined fits of 
radial q-scans in the four main directions [111], [111], [001] 
and [110] that characterize the magnetic cubic anisotropy on 
the ring (sphere). However, Eq. ([T| can be used to obtain 
the entire intensity of the fluctuations in all directions on the 
sphere. In order to verify how well our fit describes the in- 
tensity on the entire sphere we have performed a calculation 
of these intensities based on Eq. ([T| using the corresponding 
parameters obtained in the fit (see Fig.|4]in main text for the 
fit parameters). As illustrated in Fig. [8]for four representative 
temperatures the calculated intensity of the magnetic fluctu- 
ations is in excellent agreement with the experiments. The 
only temperature where major disagreement is observed is 
T = 29.1 Kjust above T c (cf. calculation I in Fig. [8]). As re- 
ported in the main manuscript between T c and 29.2 K an ad- 
ditional Gaussian shaped intensity component was observed 
at the position of the helical satellites below T c indicating 
that remanent droplets of helimagnetic order survive above 
T c . Adding these Gaussian components according to the fit 
parameters obtained via the fit of radial q (assuming isotropic 
Gaussian linewidth) for the calculation for T = 29.1 K results 




FIG. 6: The angular sections used for the integration of the magnetic 
intensities observed on the magnetic satellite peaks associated with 
the helimagnetic order in MnSi (red solid lines) and magnetic fluc- 
tuations observed on a ring (sphere) above T c — 29 K (white broken 
lines) are illustrated. Here (a) and (b) show examples of the inten- 
sity in the helimagnetic phase below T c and the fluctuations above 
T c , respectively. The black dashed line indicates a direction along 
which magnetic Bragg peaks are observed due to double scattering, 
and which has therefore not been included in the integration (see text 
for details). 
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) «(A ) 

FIG. 7: Radial q-scans through the helimagnetic Bragg peaks below 
T c = 29 K are shown for all measured temperatures. The inset in 
each panel shows the corresponding direction of the scan in recip- 
rocal space. The q-scans were extracted from the two-dimensional 
detector images in Fig.[T]in the main text by performing radial bins 
with an azimuthal width corresponding to the experimental resolu- 
tion. The solid lines are fits to Gaussian lineshapes (see text). The 
grey shaded region shown in each panel denotes the experimental 
resolution. 



in excellent agreement for this temperature as well (cf. calcu- 
lation II in Fig. [8}. 



II. GINZBURG-LANDAU THEORY OF CHIRAL 
MAGNETS 

The classical Ginzburg-Landau theory of a chiral magnet 
is given by the free energy functional T = J d 3 x f with the 
density / = f + / cub where 120 

1 ii 
fo = 2<£(r - JV 2 )0 + D0(V x </>) + -(</> 2 ) 2 - 7 0H. 

(6) 

We choose dimensionless units for the three component order 
parameter field <f>, which we normalize such that the coupling 
to the magnetic field H is just given by the magnetization 
density 7 = /is/f.u. corresponding to a single Bohr magne- 
ton per formula unit that is f.u. = 24.018 A 3 for MnSi. The 
absence of inversion symmetry allows for the Dzyaloshinskii- 
Moriya (DM) interaction D that couples internal magnetic 
space and real space and is proportional to the strength of the 
spin-orbit coupling Aso D ~ C(Aso)> which is parametri- 
cally small in MnSi. The terms of second and higher order 
in spin-orbit coupling are contained in / cu b that, in particular, 
break the rotation symmetry of fo present at H = due to 
cubic anisotropics. We restrict ourselves in the following to a 
single, representative term 

/cub = ^ {(d X <f> X ) 2 + {dy<t>yf + {d^f) + ... (7) 

where the coupling constant J cu b ~ 0(A| ). 



A. Mean-field theory 

The benchmark in the following will be the mean-field ap- 
proximation that we review here. The Ansatz for a single con- 
ical helix reads 

0hel(V) = 0o0o + ^hcie-e^ 1 - + * hcl e+e^^ r , (8) 

where 0o is the homogeneous magnetization, ^^el is the com- 
plex amplitude of the helical order characterized by the pitch 
vector Q. The vectors are given by = (ei ± ie^j^Jl with 
the normalized dreibein e\ x 62 = Q where Q = Q/Q. Eval- 
uating the energy density / with this Ansatz one obtains the 
mean-field potential V = Vo + V cu b where 

Vo=^4>o + (r + JQ 2 - 2DQ)\y hcl \ 2 (9) 

determines the strength of the amplitudes. The second con- 
tribution V cu b reflects the competition between the magnetic 
field and the cubic anisotropics to orient the pitch vector Q 
and the homogeneous magnetization (f>o, 

V cub = J<PoH(l - $ H) + ||* hcl | 2 ^(Q x O ) 2 (10) 

+ %^l*h Cl | 2 Q 2 (i-(^ + g* + ^)). 

In the paramagnetic phase | ^hci | = 0, one has (f>o = H and 
minimization <9^ V = results in the magnetic equation of 
state, 

r<Po + ^l = lH. (11) 

On the other hand, at zero field H = the homogeneous 
magnetization vanishes <fi Q = 0. From Eq. ( fTO] ) follows that in 
the helimagnetic ordered phase l^hcil > 0, even a tiny cubic 
anisotropy J cu b is sufficient to orient the pitch vector along a 
cubic symmetry direction. We assume J cu b < as this prefers 
the (111) direction as observed in MnSi [2, 3|. Minimization 
with respect to the pitch length then yields Q — D/J c h where 
JcS = J + </cub/3, and the effective potential for H = then 
reduces to 

V = <5|* hc i| 2 + ^|* hcl | 4 , (12) 
where we introduced the helimagnetic tuning parameter 

5 = (r-J cfi Q 2 )^r-JQ 2 . (13) 

In the last equation we neglected small corrections of rela- 
tive order \J cu b\/J "C 1- Helimagnetic order develops in the 
form of a second-order mean-field transition if the tuning pa- 
rameter S becomes negative at H = 0, i.e., if the parameter 
r is reduced below the DM energy density r < JQ 2 . Min- 
imization of the effective potential then yields for the ampli- 
tude l^hcil 2 = —35/u. If the magnetic field H is increased 
within the helimagnetically ordered phase and points away 
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T=29.1K T=29.3K T = 29.6 K T = 29.9 K 




Q x (A *) q x (A 1 ) 



FIG. 8: A comparison between the calculated (two upper rows) and measured intensities (bottom row) of the magnetic fluctuations above 
T c = 29 K. Here the calculated intensities where obtained by using Eq. {T} in the main manuscript that describes magnetic fluctuations arising 
at the phase transition between the helimagnetic and the paramagnetic state in a Dzyaloshinskii-Moriya helimagnet. For the calculation the fit 
parameters (s. Fig. 3 in main text) obtained by fits of radial q-scans in the four main directions [111], [111], [001] and [110] (s. Fig. 4(b) in 
main text) have been employed. For calculation II additional Gaussian profiles have been added at the position of the helical Bragg satellites 
in the helimagnetic phase for T = 29. 1 K (see text for details). 



from (111), the competition between the cubic anisotropy and 
magnetic energy in V cu b results in a reorientation of the pitch 
at a critical field iJ cl whose precise value depends on the ori- 
entation of the applied field, H, with respect to the crystallo- 
graphic (111) direction. 



B. Fluctuation spectrum and crossover scales 

The fluctuations around the mean-field solu- 
tion are described by the susceptibility tensor 
Xoij( r i r ') = /{5(pi(r)84)j(r')) with the action 
S[4>] = J d 3 xf /(ksT), which reads explicitly 

1 



Xo 



«(r,rO = 



k B T 
(0(r)0(r)^ 



(r - JVi)Sij - 2De ijn V 7 



(14) 



■24>i{r)4>j{r)) 



8{r 



0+XcuW( r ' r ') 



where 4>{r) is the mean-field order parameter and x C ub 1S 
the contribution arising from the cubic anisotropies / cu b of 
Eq. 0. 

In the paramagnetic phase, |\&hci| = 0, the Fourier trans- 
form of the susceptibility, Xq \k>k') = Xq (k)§k,-k'> 1S 
diagonal in momentum space. As shown by Grigoriev et al 
[|4], inverting this susceptibility, %q 1 (fc), for magnetic field 



H = and taking into account the cubic anisotropy J cu b 
close to the critical singularity yields the form of the scatter- 
ing cross section of Eq. ([l} in the main text, that were used 
to analyze the neutron scattering data with the identification 

"cub = I -/cub 1/(2 J). 

It is instructive to consider the mean-field fluctuation spec- 
trum wo.fc that follows from the eigenvalue equation 

fe B T^Xo]#,fc'K- (-*;') = "o, k Vi(k) (15) 



with eigenvectors Vj(k) as a function of momentum k. This 
spectrum changes qualitatively as the transition is approached 
from high temperatures, j^hell = 0, and it exhibits a series 
of crossovers as illustrated in Fig. [TJa) in the main text. Far 
above the transition, the spectrum is approximately that of a 
ferromagnet above its Curie temperature cjn.fc ~ r + Jk 2 . 
As the temperature is lowered and r reaches the order of the 
Dzyaloshinskii-Moriya energy density 



£dm — JQ 2 



(16) 



fluctuations start to become soft on a sphere in momentum 
space with radius Q and the magnetic correlations develop an 
oscillating component. The corresponding spectrum has three 
branches oj^ k < k < us^ k . The low-energy part has the 



15 



Brazovskii-form 



J 0,k 



= 5 + J(\k\-Qf 



where we neglected corrections arising from cubic 
anisotropies for simplicity. Moreover, we used the equation 
( ' of state IvPhcil 2 = — 3<5/u with the tuning parameter S of 



where 6 is the helimagnetic tuning parameter of Eq. ( fT3) . 
The dispersion of the remaining two fluctuation branches are 



Q,k 



6 + JQ 2 



Jk 2 and u 2 k 



5 + J{\k\+Q) 2 . Close 
to criticality as 5 — > + , these other modes have still a gap 
on order the DM interaction energy density £dm = JQ 2 - 
The correction to the critical spectrum ( fT7| i arising from cubic 
anisotropies of Eq. Q finally becomes important for values of 
the tuning parameter on the order of the cubic energy density, 
S ~ e C ub, with 



£"cub 



a cub £ DM = \ Jcub\Q 2 



(18) 



The energy scales that determine the crossover in the fluctu- 
ation spectrum are converted into the length scales (dm and 
( cu b of Fig.[TJa) in the main text via £dm = ^(27t/£dm) 2 and 
e cu b = J(27r/£ cub ) 2 . 



C. Mean-field magnetic susceptibility 

For zero momenta the susceptibility matrix (jT4j reduces to 
the thermodynamic dimensionless magnetic susceptibility via 

the relation Xmag,ij = &Xo,«(0,0). 

In the paramagnetic phase, | ^>hd | = 0, at zero field, H = 0, 
one obtains Xmag ij = <%(7 2 / i o)/f- In the helimagnetically 
ordered phase, j^heil > 0, on the other hand, the mean-field 
order parameter 4>(r) carries a finite momentum Q and, as a 
result, the susceptibtility is non-diagonal in momentum space. 
Using the Ansatz ([H} we obtain 



Xo,zj(k, k ) 

" 2 



Xolji k ) S 0,k- 



<5_ 



(19) 



3 



2Q.k+k' 



*2 + 



her 



6j "2Q,fc+fc' 



Q,k+k' 



The part that is diagonal in momenta reads 
1 



xoi,( fc 



k R T 



(r + Jk 2 )S, K 



2D *£{j n zk ri 



Eq. ( fT3| l that is negative in the helimagnetically ordered 
phase, 5 < 0. Close to the transition 6 — > 0~ the magnetic 
susceptibility becomes isotropic and smoothly connects to the 
value in the paramagnetic phase. Deep in the ordered phase, 
S <C — JQ 2 , the susceptibility is anisotropic. In particular, the 
longitudinal susceptibility, Xmag = Xmag,ijHiHj, depends 
on the angle between the magnetic field and the pitch vector, 
HQ, 



mf _ 7^o (1 - S/(JQ 2 )) + (HQ) 2 (-S/(JQ 2 )) 

Amag ~~ T ^i9 i nsr// t^->9\ ' 



JQ 2 



1 - 25 /{JQ 2 ) 



The macroscopic magnetic susceptibility (xmag) averages 
over helimagnetic domains that might possess pitches in dif- 
ferent (111) directions. For field cooling (FC) with H in a 
(111) direction, we might expect a single helimagnetic do- 

main with (HQ) 2 = 1 so that (xZg)^ = tWW) 
attains the maximal possible value independent of 5. This is 
also the value of the longitudinal susceptibility one expects in 
the conical phase at a finite field H > H c \. For zero field 
cooling (ZFC), on the other hand, we can assume an equal 
distribution of domains so that ((HQ) 2 ) = 1/3 independent 
of the magnetic field orientation 



(x MF )\ 

VAmag/ I 



ZFC — 



7V0 1 - p/(JQ 2 ) 
JQ 2 1 - 28 /(JQ 2 ) ' 



(23) 



As the helimagnetic phase is entered this susceptibility drops 
from the mean-field value at criticality "f 2 Ho/(JQ 2 ) and sat- 
urates for S <C — JQ 2 to a value that is reduced by a factor of 
2/3. 



■ 2<f)Qi4>Qj) + -|V'hci| 2 (2% - QiQj] 



The magnetic susceptibility is obtained from the inverse 
of the generalized matrix Xoy(k>&') taken at zero mo- 
menta. It is however important to note that Xo,ij(0, 0) 7^ 
(Xoij(®i 0)) _1 because the order parameter carries momen- 
tum. After inverting the susceptibility matrix, we get for 
the dimensionless magnetic susceptibility of a single helimag- 
netic domain Xma g ,ij = aT^t Xo,ij (0: 0) at zero magnetic field 



HI. BRAZOVSKII THEORY FOR CHIRAL MAGNETS 

Close to the critical temperature, chiral magnetic fluctua- 
tions become very important and have to be treated in a self- 
consistent manner. Corrections arising from the interaction of 
(20) modes with the dispersion of Eq. (TJ\ are singular as these 
fluctuations have a quasi one-dimensional character. As the 
XcubijW -momentum dependence is peaked at a finite momentum along 
the radial direction, the density of states 



v(e) 



(2tt) 3 o,kJ 2ir 2 ^/jQ 2 Ve-5 

(24) 



H = 0, i.e. 



= 



MF _ 7 Mo 

A mag, i j 



JQ 2 



l-S/(JQ 2 ) 
1 - 25 /( JQ 2 ) 



iQi) 



(21) 



possesses a one-dimensional singularity as e approaches the 
tuning parameter, e — > S. Brazovskii [6| argued that the in- 
teraction correction due to fluctuations with such a singular 
density of states changes the second-order mean-field transi- 
tion into a fluctuation-induced first-order transition. 

In the following, we discuss in detail the application of 
Brazovskii theory to chiral magnets. We start from the two- 
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particle irreducible (2PI) effective action |7-9| 

m,x]= (25) 

= S[<P\ + ^Trlogx- 1 + ^Tr (x^ 1 - x" 1 ) X + X ], 

which is a functional of the field configuration <fr and the prop- 
agator x- The action S[<p] = J-[<p]/ (k B T) is given by the free 
energy functional of the chiral magnet defined in the previous 
section. Here the trace, Tr, should be taken in the functional 
sense. The bare susceptibility was already defined in Eq. (14) . 
Finally, T2 is the sum of all 2PI vacuum graphs evaluated with 
the renormalized propagator x, see Fig. [9] The effective ac- 
tion has to be minimized such that 



ar[0,x] 

5(f>(r) 



= 0, 



ST[<P,X} 
8x(r,r') 



= 0, 



(26) 



which determines the observable field configuration <fi and the 
susceptibility x- 



A. Hartree-Fock-Brazovskii (HFB) approximation for the 
effective potential 

The following approximation consists of two steps. First, 
we will limit ourselves for the 2PI vacuum graphs to the low- 
est order diagram in Fig. |9ja) corresponding to the Hartree- 
Fock approximation, 



r. 



A\k B T 



dr 



(trx(r,r)) 2 + 2tr X 2 (r,r) 



(27) 



The second approximation concerns the form of the suscep- 
tibility matrix. It was argued by Brazovskii [6 | that close to 
the transition the most singular contribution of the fluctuation 
correction is attributed to the part of the fluctuation propaga- 
tor that is diagonal in momentum. We thus make the following 
Ansatz for the susceptibility 

X^(k,k') =x^(fc)<Wfe' +Xojj( fe > fc/ )( 1 - <Wfe') 

(28) 

While the off -diagonal part coincides with the one of the bare 
susceptibility ( fT9| ), we parametrize the diagonal part as 



Xij 



(*0 



k R T 



(JQ 2 + Jk 2 )d l3 
■FAj_ (Sij - hifij) + Aynifij 



(29) 



with the two variational parameters Aj_ and A|| , and the unit 
vector h. For our purposes it is sufficient to assume that for 
Aj_ A|| there is only a single preferred orientation given by 
n. At any finite field H in the paramagnetic phase j^hell = 
or in the helimagnetically ordered phase, | "Jhci | > 0, for fields 
H > H c i, the unit vector n can be identified with the mag- 
netic field orientation, n — H. For l^hcil > and H = 0, 
on the other hand, the orientation is determined by the pitch 
orientation of the helimagnetic domain so that n = Q. How- 
ever, for H H c i within the ordered phase where cubic 



anisotropies compete with the magnetic energy to orient the 
pitch, see Eq. ( (TO) , two orientations might be present and the 
parametrization of Eq. |29) might be insufficient. 

Confining ourselves to the critical regime and anticipating 
that the crossover close to criticality associated with the en- 
ergy scale e cu b of Eq. ([18) is basically preempted by the first 
order transition, we neglect in the following the corrections 
from cubic anisotropies to Eq. (129) , Furthermore, we apply 
the Brazovskii approximation J6J 



Xo 



y (fe, fc ) 



Xoi)( k ) s o. 



k+k' 



x^(fc,fe') 



(30) 



in the evaluation of the terms in ((25), i.e., we neglect the com- 
ponents off-diagonal in momentum space which are present 
in the ordered phase. Eq. ((30) together with Eq. ((27) yields 
the self-consistent Hartree-Fock-Brazovskii (HFB) approxi- 
mation for the effective potential. 

For the evaluation of Eq. ( |25) , it is convenient to define the 
function 



£>(A_l,A||) = k B T 



dk 



\ogx-\k), (31) 



with the susceptibility of Eq. ((29). In the following, we 
will also need derivatives of this function which we de- 
note by a corresponding subscript, e.g., D_l(A_l, A|| ) = 
9a ± D(A±, Ay). The function D can be separated into a part 
containing the leading singularity and a part, that is sublead- 



ing close to criticality, D 
is given by 



-D su b- The singular part 



D sing (A ± ,A|, 



Q 3 k B T / Aj_ + Am /Aj_ — A 



V2ir V J Q 2 
where the auxiliary function y reads explicitly 



■An 



(32) 



y(a) = n/TT^ 



ds 1 



1 + s 2 



arctan 




and obeys y(0) = 1 and y'(0) = g. The subleading part has 
the property that its first derivatives have a well-defined limit 

as A_l,|| 0. 

Using the mean-field Ansatz |8) for the field configuration, 
the effective action V then reduces to the effective potential 

r = vv cS /(k B T) with 

V eff (A X , A|| , 0o, Vhcl, ^hel) = Wo, V>hel, V'hcl) + ~-D(A_L, A||) 



+ 1 -(s+^l + 4\^\ 2 )-A A _)D ± 

" / ■ " \2 , of 1 r>2 , n 2 



(34) 



where Vo is the mean-field potential given in Eq. |9). We as- 
sume that the part of the mean-field potential V cu b of Eq. ( fLO) 
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determining the orientation of the pitch Q and the orientation 
of the homogeneous field (j>o is already minimized and that in 
the ordered phase it yields (f>o = Q for finite </>n, which might 
exclude the regime close to the pitch-flop transition that we do 
not consider in the following for simplicity. 

The effective potential V G ff now is to be minimized with 
respect to Aj_ and Am, 0o and the complex amplitude V'hei- 
Minimization with respect to the susceptibility parameters, 
3a ± V e ff = and 9a » Veff = yield the two equations 



A ± =5 +3! , 



+ 4|tM 2 ) + ^(2Z? ± + /J||) (35) 

Ay = 8 + |(3$ + 2|^ hcl | 2 ) + | (D ± + 3D||) (36) 

Minimization with respect to the amplitudes, d<f, V c s = and 
d,p* j V e fi = gives the equation of states 



(05 + 2|tM 

2 i o|j, |2 



3! 

u 

+■ 3j(0S + 2 l^hci 



2D ± 



-3D||)j0 o = 7 H (37) 
= 0. (38) 



Taking the effective potential V c ff at its minimum yields the 
free energy density in the HFB approximation from which 
thermodynamic quantities can be evaluated. A numerical cal- 
culation of the specific heat is shown in Fig. [5] in the main 
text. One finds that the singular fluctuation corrections renor- 
malize the effective potential strongly so that the second-order 
mean-field transition is converted into a fluctuation-induced 
first-order transition. 



(a) 




(b) 



FIG. 9: Two-particle irreducible diagrams of first (a) and second (b) 
order in the interaction. The interaction is represented by the dot and 
the line corresponds to a fluctuation propagator. We approximate T2 
of Eq. {15) with the lowest order diagram (a). 



B. Magnetic susceptibility at zero field H=0 



With the help of the susceptibility tensor d28l we can de 



rive the magnetic susceptibility, Xmag.ij = \^rXij (0>0). 
following the derivation of the previous section. As before, 
we limit ourselves to the magnetic susceptibility in zero field, 
H = 0. In the paramagnetic phase, T > T c , the magnetic sus- 
ceptibility is isotropic Xmag,ij = 8ijj 2 /iq/(JQ 2 + A), with 



A = A || = Aj_. In the ordered phase, T < T c , we obtain for 
the magnetic susceptibility of a single domain 



Xmag,ij — 

7 2 M0 / QiQj 

JQ 2 Vl + A|,/(JQ 2 ) 



1 + A ± /(JQ 2 ) 
1 + 2A ± /(JQ : 



(39) 



where we used that Aj^ = (w/3)|^hci| 2 which follows from 
combining the equation of state ( (38] l and Eq.([35]>. Averaging 
over domains as before, we thus obtain for the longitudinal 
susceptibility in the case of zero field cooling (ZFC) 

(40) 

for T > T r 





1 + for T<Tc 

A||/£ D M 1 + 2A_l/£ DM / 

where £dm = JQ 2 , and A^, A|| are functions of the tuning 
parameter 8 that have to be determined by minimizing the ef- 
fective potential. Far-away from the transition the mean-field 
behavior is recovered but close to the transition the Brazovskii 
renormalization leads to qualitative modifications as discussed 
in the main text. Fig.|4|d) in the main text shows a comparison 
with the experimental data on MnSi. 



C. Brazovskii scaling limit in the critical regime 

In the presence of fluctuation corrections an additional en- 
ergy scale £Qi ~ (ufcsTQ 2 /v / J) 2 ^ 3 emerges. As the tran- 
sition is approached from high temperatures at H = 0, the 
fluctuation corrections become non-perturbative if the tuning 
parameter 8 reaches the order of 8 ~ egi, which corresponds 
to the Ginzburg criterion for chiral magnets (provided that 

£Gi < EDM)- 

It turns out that in the critical regime the effective poten- 
tial acquires a scaling form if the Brazovskii energy is much 
smaller than the DM energy density, Eq; <C £dm- This scal- 
ing limit is obtained when the function D is approximated by 
its singular part D S i ng only. For example, in the paramagnetic 
phase, iV'hcil = 0, at zero field, H = (j>o = 0, the equations of 
state are automatically satisfied, and Eqs. ( [35] l and ( [36] l reduce 
to a single equation as Aj^ = Am = A that has a simple form 
in the scaling limit given by 



3/2 

£ Gi_ 

/A' 



(41) 



This equation has the Brazovskii-form; the 1/vA depen- 
dence of the fluctuation correction reflects the singular density 
of states ( |24] i. Eq. ( |4T] i defines the Ginzburg energy density 
including prefactors to be 



£Gi 



5 uk B TQ 3 

7W 



2/3 



36tt 



With the identification A = Jn 2 , 8 
Eq. ( |4Tj ) reproduces Eq. ([3]) in the main text. 



Jk^p andf 



Gi 



(42) 
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D. Corrections to the Hartree-Fock-Brazovskii approximation 

The HFB approximation consists of (i) approximating the 
2PI vacuum graphs only by the lowest order one, Eq. $2J\ , and 
(ii) neglecting the off-diagonal components of the susceptibil- 
ity matrix, Eq. ( f3Q] >. In the following, we discuss the validity 
of these two approximations separately. 

The correction to the Hartree-Fock approximation of T 2 
corresponds to the next-to-leading diagram shown in Fig.[9jb). 
In the scaling limit, we find that this diagram results in a sin- 
gular correction to the right hand side of Eq. ( |4T| > of the order 



(43) 



In order to estimate the regime where this correction is negli- 
gible we follow Ref. [6| and demand that <5A/A -C 1. Close 
to the fluctuation-induced first-order transition at H = 0, 
the parameters are determined by the Ginzburg scale so that 
A ~ g Gi and this condition becomes 



SA 

aT 



£Gi 
EDM 



< 1 



(44) 



The approximation (i) is thus self-consistent as long as the 
Ginzburg scale is much smaller than the DM energy scale. In 
the case of MnSi, we find kgi/ k dm — ^dm/£,Gi ~ 0.5, so 
that quantitative corrections due to Fig.[9|b) can be sizeable. 

The susceptibility matrix, Eq. ( [30) 1, is non-diagonal in mo- 
mentum space as the helical order parameter carries momen- 
tum. The approximation (ii) neglects the off-diagonal com- 
ponents of the susceptibility within the ordered phase. This 
approximation will account for the quasi- Id renormalization 
arising from the Brazovskii-spectrum close to criticality but it 
neglects corrections arising from the Goldstone mode within 
the helimagnetically ordered phase. These latter corrections 
can become logarithmically large in the isotropic limit if 
H = and cubic anisotropics are absent. The dispersion of 
this Goldstone mode is then anomalously soft ifTUl ITTI result- 
ing in the absence of true long-range order as a consequence 
of the so-called Landau-Peierls instability familiar from smec- 
tic liquid crystals fl2l . Importantly, the cubic anisotropics in 
the present case will regularize the logarithmic singularities 
so that true long-range order sets in at a finite critical tem- 
perature. We assume that the cubic anisotropics are small 
enough, e cu b <C £Gi> so that they can be neglected in the 
Brazovskii renormalization of the effective potential but are 
sufficiently large in regularizing the logarithmic singularities 
within the ordered phase. The study of the residual logarith- 
mic terms and the modifications of thermodynamic anomalies 
in the limit of vanishing cubic anisotropics will be addressed 
in a separate publication. 



IV. ESTIMATE OF PARAMETERS AND 
CHARACTERISTIC SCALES 

We discuss in the following the interpretation of thermody- 
namics of MnSi close to the critical temperature in terms of 



TABLE III: Estimate of microscopic parameters of the Ginzburg- 
Landau theory of Eqs. {6} and {7]l close to the critical temperature, 
see text 



Q R 


i 0.039/A 


J P 


d 2.8 meV/A 


|^cub| s 


i 0.13meV/A 


u R 


t 0.32 meV/A 3 



TABLE IV: Estimate of characteristic energy densities, see text. 



units of 


lCT 3 meV7A 3 


EDM 


EDM 


4.3 


1 


EGi 


1.0 


0.23 


^cub 


0.1 


0.023 



the HFB approximation for chiral magnets. In particular, we 
discuss the values of the microscopic parameters listed in Ta- 
ble III of the Ginzburg-Landau theory of Eqs. (|6]l and |7]) that 
have been used for the comparison. 

The length of the pitch Q ps 0.039/A close to the transition 
temperature was determined with the help of our SANS data, 
see Fig.|4jb) in the main text. From the fit of the magnetic sus- 
ceptibility to Eq. |4} of the main text we obtained \o ~ 0.27 
that can be identified with \o = Y 2 Mo/(JQ 2 )> from which 
follows for the stiffness J « 2.8meWA. This determines the 



DM energy density ebm ~ 4.3 x 10 



-3, 



zV/k 3 . From the 



fit to the magnetic intensity we extracted the cubic anisotropy 



parameter a c 



I Jcub|/(2J). It assumes a maximal value 



ib 

ib 



0.023 at the transition from which we can estimate 
■i 0.13meV7A. For the cubic energy scale we obtain 



I Jo 

£ C ub = Ja 2 cuh Q 2 ps 0.1 x 10- 3 meF/A 3 . The value for the 
interaction u is more difficult to determine. One possible way 
is to use the magnetic equation of state of Eq.([TTJ. Compar- 
ing with magnetization measurements |5| on MnSi at small 
fields along (100) at a temperature T — 32 K we obtain for 
the interaction constant UAnot ~ 0.13meWA 3 . However, 
experimentally the magnetic equation of state has not the sim- 
ple mean-field form ( fTT| i probably due to strong Brazovskii 
renormalizations and the error for the extracted value of u is 
probably large. For this reason, we used instead the value that 
follows from the Brazovskii fit to the correlation length that 
yielded a Ginzburg length Kq{ pa 0.019/A. This translates to 



a Ginzburg energy density eci 



Tk 2 
Gi 



1.0x10" 



2V/A 3 

and, in turn via Eq. ( |42) , to the value for the interaction 
u ps 0.32mey/A 3 . 

Finally, the temperature dependence of the tuning param- 
eter was extracted via the Curie-Weiss limit of the magnetic 
susceptibility. Experimentally, one finds sufficiently far above 
the critical temperature that the dimensionless magnetic sus- 
ceptibility obeys Curie- Weiss behavior Xmag = T X /(T — T c ) 
with a Curie temperature 



T c ps 28.2^ and T x 



0.62X. 



(45) 



The values were extracted from the asymptotic behavior of the 
fit to the susceptibility in Fig. |5]l(a) in the main text. The tem- 
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perature dependence of the tuning parameter S oc T — Tmf or, ered. This fixed all microscopic parameters and the resulting 
equivalently, r of Eq. |6]) was determined by adjusting the sus- HFB approximation yielded the thermodynamics displayed in 
ceptibility so that the observed Curie- Weiss behavior is recov- Fig. Q in the main text. 
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